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The phenomenology of the five independent sets of 3 x 3 quark mass matrices with five texture 
zeros is 

carried through in full detail, including predictions for the CP violation asymmetries. Our study 
is done without 

any approximation, first analytical and then numerical. 



I. INTRODUCTION 

Although the gauge boson sector of the Standard 
Model (SM) with the SU{i)c ® SU{2)l ® U{1)y local 
symmetry 

has been very successful its Yukawa sector is still 
poorly understood. Questions related with this sector as 

for example the total number of families in nature, 
the hierarchy of the charged fermion mass spectrum, the 
smallness of 

neutrino masses, the quark mixing angles, the neutrino 
oscillations and the origin of the CP violation, remain 
until today 

as open problems in theoretical particle physics. 

The mechanisms for fermion mass generation and fla- 
vor mixing can be classified into four different categories: 
(i) radiative mechanisms (which includes the so-called 
Froggat and Nielsen mechanism Q); (ii) texture zeros in 
the mass matrices which may predict self-consistent and 
experimentally favored relations between fermion masses 
and flavor mixing parameters [3,il; (iii) family symme- 
tries, discrete @ and continuous, global and local gauge 
symmetries [7|], and (iv) see-saw mechanisms for electri- 
cally neutral Q and charged particles @, related to a 
natural interpretation of the smallness of some fermion 
masses. 

In the SM and after the local gauge symmetry has been 
spontaneously broken, the quark mass terms are given by 



Cm = UolMuUor + DoLMdDoR + h.. 



(1.1) 



D, 



OL 



_where Uql = (uo,co,to)L, 
(do. So, 6o)l, Uqr ~ (uo,co,io)fi, Dqr = (do, sq, 69)^, 
(the upper T stands for transpose, and the down zero 
stands for weak basis states). The matrices M„ and Md 
in p. II) are in general 3x3 complex mass matrices. In 
the most general case they contain 36 free parameters. 
In the context of the SM, such a large number of 
parameters can be drastically cut by making use of the 
polar theorem of matrix algebra, by which, one can 
always decompose a general matrix as the product of an 
hermitian times a unitary matrix. Since in the context 
of the SM the unitary matrix can be absorbed in a 



redefinition of the right handed quark components, this 
immediately brings down the number of free parameters 
from 36 to 18 (the other eighteen parameters can be 
hidden in the right-handed quark components in the 
context of the SM and some of its extensions, but not in 
its left-right symmetric extensions). 

So, as far as the SM is concerned, we may treat without 
loss of generality M„ and as two hermitian quark 
mass matrices, with 18 real parameters in total, out of 
which six are phases. Since five of those phases can be 
absorbed in a redefinition of the quark fields [l3|, the 
total number of free parameters we may play with in M„ 
and Md are 12 real parameters and one phase; this last 
one used to explain the CP violation phenomena. 

But in the context of the SM it is always possible to 
implement the so called weak basis (WB) transformation, 
which leaves the two 3x3 quark mass matrices Hermitian, 
and do not alter the physics implicit in the weak currents. 
Such a WB transformation is a unitary transformation 
acting simultaneously in the up and down quark mass 
matrices \1m- That is 



M. 



UMuU\ 



(1.2) 



where U is an arbitrary unitary matrix. We say then 
that the two representations {Mu,Md) and {M^,M^) 
are equivalent in the sense that they are related to the 
same Vck m mixing matrix. This kind of transformation 
plays an important role in the study of the so-called flavor 
problem. 

In the last paper of Ref. [ll| it was shown that, re- 
lated to the mass hierarchy m„ < rric < mt and < 
TOs < rrib, it is always possible to perform a weak ba- 
sis transformation in the hermitian quark mass matrices 
such that 



(A/«)i3 = {M^)3i - 0; 

or equivalently; (1.3) 

(Mf )n = (Af«)n = (M«)i3 - (Mj^Oai - 0. 

According to this, it is always possible to have hermi- 
tian quark mass matrices with three texture zeros which 
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do not have any physical imphcation. With three texture 
zeros the number of free parameters in and re- 
duces from twelve to nine real plus one phase, just enough 
to fit the measured values for the six quark masses, the 
three mixing angles, and the CP violation phenomena. 
Any extra texture zero can only be a physical assump- 
tion and should imply a relationship between the quark 
masses and the parameters of the mixing matrix. But the 
maximum number of such texture zeros consistent with 
the absence of a zero mass eigenvalue and a non degen- 
erate spectrum is just six, with three in the up and three 
in the down quark sectors respectively. 

In what follows we are going to present analytic and 
numeric results for the set of SM quark mass matrices 
containing five texture zeros, taking special care to ac- 
commodate the latest experimental data available [l^ . 
including the CP violation phenomena. 

This paper is organized as follows: in Sec. [IT] some fea- 
tures of the SM mixing matrix are presented, in Sec. IIIII 
we introduce the five independent sets of 3 x 3 quark mass 
matrices with five texture zeros, in Sec. IIVI we present 
our analytic and numerical analysis. A discussion of the 
results is done in Sec. |V] and the conclusions are pre- 
sented in Sec. |Vll Two appendices are written at the 
end, the first one presenting the experimental measured 
values used in the main text, and a second one is a pure 
mathematical appendix, dealing with the analytic diag- 
onalization of real 3x3 orthogonal matrices. 



II. THE SM MIXING MATRIX 

In the SM and for the six flavor case, the Baryon 
charged weak current is given by 



Vc 



KM 



Vud Vus Vub 
Vcd Vcs Vcb 

Vtd Vts Vtb 



(2.2) 



The unitary of the CKM mixing matrix leads to rela- 
tions among the rows and columns of VckMi in particular 
we have for the columns: 



Vudv:, + Vcdv:, + Vtdv:, = o, (2.3a) 
Vusv:^ + ycsv:^ + ytsv:^ = 0, (2.3b) 

KdK*6 + KdK; + = 0. (2.3c) 



Each of these three relations requires the sum of three 
complex quantities to vanish and so can be geometrical 
represented in the complex plane as a triangle. These 
are the unitary triangles though the term "unitary 
triangle is usually reserved for the relation (|2.3cp only. 

The three angles of the unitary triangle represented by 
(j2.3cp , which are physical quantities and can be indepen- 
dently measured by CP asymmetries in B decays, are 
defined as follows llj: 



arg 



/3 = arg 



7 = arg 



th 



VcdVX 



VtdV,l_ 

V^udVub 

VcdV;, 



(2.4a) 
(2.4b) 
(2.4c) 



(2.1) 



where Vckm = UuUl is the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix, with [/„ and Ud the 
unitary matrices which diagonalize the Hermitian M^MI^ 
and MdM^ square mass matrices respectively, and Ul = 
{u,c,t)]^ and = {d,s,b)J^ stand for the quark field 
mass eigenstates. 

Vckm is a 3 x 3 unitary matrix, its form is not unique, 
but the permutation freedom between the three genera- 
tions can be removed by ordering the families such that 
(ui,W2,U3) — {u,c,t) and {di,d2,d3) — > {d,s,b). The 
complex elements of Vckm are thus commonly written 
as 



The experimental findings at the B factories, fitted to 
close the triangle, are [13, [3 



(a, /3, 7){;*„ = (95.9 



-1+2.2 
-5.6' 



21.8 ±2. 



67.2 



+4.4- 
-4.6, 



, (2.5) 



with an accuracy in the measurement of sin 2/? no less 
than 20% O. 

The Cabbibo-Kobayashi-Maskawa Vckm matrix, can 
be parameterized by three mixing angles 9, 



< 



j, i, J = 1, 2, 3, the angles between the i and j fami- 
lies, and only one CP violating phase [l3|- Of the many 
possible paramcterizations, the standard choice is [l2l.[T7| 



C12C13 S12C13 sise"*"^ 

Vckm = | -S12C23 ~ 212523^136*'' C12C23 - si2S23Si3e"' S23C13 | , (2-6) 

S12S23 - Ci2C23Sl3e'* -C12S23 - Si2C23Sl3e*'' C23C13 
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but the most important fact related with this matrix 
is that most of its entries have been measured with high 



J 



accuracy, with the fohowing 
bounds [Hill: 



0.970 < \Vud\ < 0.976 0.222 < |Ks| < 0.226 0.003 < |Kb| < 0.004 
yiexp) ^ I Q < iv^^i < 0.237 0.960 < \Vcs\ < 0.990 0.039 < |Kb| < 0.041 
0.008 < l^tdl < 0.009 0.038 < l^t^ I < 0.042 0.999 < |Vth| < 1.000 



r 



(2.7) 



where the experimental numbers quoted above at 95% 
C.L. are beyond the pure experimental bounds (which 
caU for example for a \Vtb\ > 0.79) because they are 
restricted to fit the unitary conditions of Vckm- The 
numbers quoted in (j2.7p are the most convenient for our 
purpose, due to the fact that we are going to confront 
these numbers with quark mass matrices which must fit 
the SM constraints. 



Pattern 



III. FIVE TEXTURE ZEROS 



As mentioned before, the maximum number of texture 
zeros in the quark mass matrices of the SM, consistent 
with a non degenerate spectrum and with the absence of 
a zero mass eigenvalue, are three in the up and three in 
the down quark sectors respectively. With this in mind, 
Harald Fritzsch proposed some time ago [3| the existence 
of parallel three texture zeros structures for the 3x3 
quark mass matrices M„ and Md of the SM, such that 
(M,)n = (Af,)22 = (M,)i3 = (M,)3i =0, q ^ u,d. 
This original Fritzsch "anzatze" is named today in the 
literature as the "parallel nearest neighbor interaction 
form". With six texture zeros there are left just six 
real parameters plus one phase to describe the six quark 
masses, the three mixing angles and the CP violation 
phase, and so, the three mixing angles of the CKM mix- 
ing matrix can be expressed as functions of the quark 
masses and of the CP violation phase. 

As it has been shown in several places, the parallel and 
non-parallel six texture zeros does not work properly and 
it has been ruled out by analytic [l^ and numeric [l^ 
studies (only a charm quark mass half of its measured 
value can rescue this " anzatze" ) . 

Later, Ramond, Roberts and Ross (RRR) gave up the 
parallelism between the structures of M„ and Md and 
found that there exist five phenomenologically allowed 
patterns of Hermitian quark mass matrices, each one of 
them with five texture zeros as listed in the table. 



II 



III 



IV 



V 



M* 




M 



W 



bd 
dd 



ad 

fld c<i 

dd 

ad 

a*d Cd 

hi 

ad 

a2 Cd bd 

&S dd 

ad 

a*^ Cd 

dd 

fld 



Cd 



bd 
dd 



Five RRR patterns of Hermitian quark mass matrices. 



A few remarks are in order: 



i) Patterns I, II, IV and V have 3 different phases, so 
2 are unphysical. Pattern III has four different phases, 
so 3 are unphysical. 



ii) Each pattern has 7 real parameters which means 
two physical predictions. That is, for each pattern 2 of 
the 3 mixing angles can be written as a function of the 
six quark masses and of the CP violation physical phase. 



iii) The five different RRR patterns commit at least 
with one (some with both) of the weak basis texture ar- 
rangements in Eqs. (jl.3|) . 

In what follows we are going to study in detail these 
independent RRR patterns, paying special attention for 
each case to the prediction for the Cabbibo angle and to 
the CP violation phenomena. 
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IV. ANALYTIC AND NUMERIC ANALYSIS 

1. Pattern I 

This pattern was studied previously in Ref. The 
complex Hermitian quark mass matrices for this pattern 
can be written as: 



(4.1) 






















V 


|6«|e-'^" 


du 








(4.2) 



As can be seen, for this pattern the mixing angles 6113 
and 023 between the third family and the first two 

ones, come only from the up quark sector. The com- 
plex phases are removed by going to a prime basis using 
the following unitary transformations: 





where t/™^ and U^J arc complex diagonal matrices as 
presented in appendix (jB|) . 

So, in the primed basis, the algebra reduces to diag- 
onalize the two real symmetric mass matrices (jB17l) for 
q = u 

and (|B7p for q = d, which are: 




Mi 



\ad\ 
Cd 
dd 



diagonalization carried through in Appendix (|B|) . 
From the former analysis we can evaluate the mixing 



matrix Vqj^j^j 



0: 



bfiT 



, where and t/, 



bf 



are as defined in (|4.3p and calculated in Appendix ([B]) 
and O^^ and O^/ arc given by (jB19|) and (|B9p respec- 
tively. 

The elements of the CKM mixing matrix can be ex- 
pressed in the context of this pattern as 



CKMH^. = {0:'')n{0'/)rnl+e^*HO:''hiO''/)rn2 

(4.4) 



+e*(0:'');3(Oy)™3 



where 0i = (a™ - a^J) and </)2 = « + Pu)- 
The use of the matrices in appendix ([B]) allow us to 
write the following analytic mixing matrix: 



CKM 



:(v^ 



du{rn^+ms){mt—rau){^c+fnu) 



^ du{md+ms)(mc+mt)(mc+mu) 
^/ — dit-\-mt(^e^'^^ \/dumdmt-\-^mcmsmu^ 
y^d^,{ma+ms){ma+mt){mt—mu) 



-y/ du{'md+ms){mt—mu){fnc+mu)- 
e^'^i y/duync^s-^-^/^d^t^u 

/ d„(md + ms)(n»c + mt)(nic + m„) 
V du+TTic 

y/d^(md+nis){ma+nit){mt—mu) 



^i02 (d„-|-Trae)Trat(d„-m„) 

d^{ma+mt){mt-mu) / 



where as anticipated, the last column is a function of 
the parameters in the up quark mass matrix only, and a 
common phase 

g»02 jj^ that column can be removed by a redefinition 
of the b quark field. 

Our approach now is the following: first we use the 
free parameter du in order to fit the experimental value 
Kb 

and then use the phase (f>i to fit | Vus \ ; when the central 
quark mass values are used we obtain for this pattern 

the result du ~ mt and 4>i ~ 1.6. The next step is to 
fine tune both du and using a random 

numerical analysis, using quark masses niu and at 
2(7 values with the other four quark masses at 

Icr [2^), aiming to get results as close as possible to 



the ones presented in F^'^^p) in (E^]). 

Finally, for the best values obtained we calculate the 
CP asymmetries a, f3 and 7 as in (|2.4p . 

Notice that a crude approximation for which du = mt, 
using further mt » mc » mu, implies from the mixing 
matrix 




a form advocated in several references [2T| . But to set 
du — mt implies bu — according to Eq. (jBlSp , and 

not only we get back to a parallel six texture zero pat- 
tern with = 623 = (something that can be seen 
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also 

from V^j^j^j, the 3x3 mixing matrix above). So, om- 
fine tunning approach, is mandatory, c?„ must be different 

from rrit and any analysis done for which expression 
(|4.5[) plays a central role, must produce dubious predic- 
tions. 



In this regard, the conclusions in Rcf. [21| . which are 
contrary to ours, may be wrong. 

The numerical subroutines used throws for pattern I 
the following numbers (mass parameters are in GeV's 
and angles in radians): 



0.0023, 



0.560, mt 172, 



rud = 0.0029, rus = 0.06, rub = 2.89, 

du = 171.721 < mt, (/-i = 1.61; 
which imply the following 3x3 mixing matrix: 



0.97428 0.22532 0.00264 \ 
0.22517 0.97349 0.04020 
0.00865 0.03934 0.99919 / 



which not only close the triangle, but are such that a 
and 7 agree with the measured value at Icr and 
l3 at 2(7. 



2. Pattern II 

Notice for this pattern that the mixing angles ^13 and 
023 between the third family and the first two 

ones come only from the down quark sector. Also, 
this pattern is obtained from the previous one by the 
replacements: 

Mj — > M^^ and — > il/j^ and d„ dd, where 
dd is the free 

parameter of this pattern. 

As in the previous pattern, we start by removing the 
complex phases from M^^ and M^^ and then use the 
results 

in appendix B in order to diagonalize the orthogonal 
mass matrices. When done, the CKM mixing matrix for 
this pattern reads: 



^CKM 



0'L^Ut^'U^''''[0^''f, where the orthogonal 
and unitary matrices are presented in appendix 

(B). The elements of the CKM mixing matrix can be 
expressed now as 



numbers in quite good agreement with the experimen- 
tal measured values. Finally, the three angles of the uni- 
tary triangle 

of the B decays CP asymmetries calculated according 
to (|2.4p for this pattern are: 



(a, /3, ^)if^ = (90.79, 16.51, 72.68), 



(VS'kmW. = {0':/)iii0f)mi+e^'^^{0'jM0fU2 
+e^^n0'/)/3(0:f'')™3, (4.6) 

where (j)i = {a'^ — a"^) and 02 — — ("d + )' which 
allows us to write the following analytic 
forms: 
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Vud 



* us 



^ ub 



V" 



II 



^ mcnib {dd - md) e"^i ^77i„ - rud) md 



\/(t71^ + Uln.) (Tflh — Tfl^) (ui^ -\- Tflc) 
y V c 1 ' ' ) V D ''"a/ V, ' s/ 


\I{tii^ + m.,^ fm?. — jtia\ (m^ 4- 


' '"s ! 


^J mcmi,md {dd + itIs) 


e"*^! mgTOu (drf + m^) 




\ ivcir- -4- Tfi^, \ dr\ (Tfih + TTiol (yriw 4- yTi^l 


\/ (nir- -\- TTi^A (uih + m^. 1 fmw + 


' ' "S J 


\/ {—dd + ?77,;,) mdmsTric 


e'"^i -v/ mb (— + mb) mu 




w 1 1 ' / ^ci \' ' ''0 ' ' ''d J \' ' '^0 1 ' ' ''S J 


w 1 1 ' ' '^U J \' ' ''0 ' ' '^d } \' ' ''0 1 


' ' ) 


e"^^ ^ {dd - md) uidmc 


^J muTrib {dd - md) rus 




^ (toc + m„) (mb - nid) (jrid + rris) 


^ {rric + uiu) dd {rub - rud) {nid 4- 


' ms) ' 




^Jmi^iTid {dd + nis) rriu 




a/ (mc + m„) (m^ + m^) (m^ + m,,) 


^ {iTLc + TTiu) dd {ruh + rris) {md 


+ ms) 


e*"^! rricmb {-dd + mb) 


^ {-dd + m^) uiumdms 




{rric + m„) (mb - m^) (mb + m^) 


y (777,c + - md) {mb + 


rris) ' 



_ ,02 / (mb - dd) TOd {dd + 

]] dd{mb - md) {md + ms)' 

/ (^fc - dd) (dd - md) ms 
]] dd{mi, + nis) {md + ms)' 

/ mb (dd - md) (dd + m^) 
y dd{mb - md) {mb + ms)' 



(4.7) 



where again the last row is a function of the parameters 
in the down quark mass matrix only, and the phase 

6**^2 can be removed by a redefinition of the field for 
the quark t. The form of the former mixing entries 

recommends to use dd to fit the experimental measured 
value Vts and then (j>i to fit Vus, which 

produces dd ~ mb and (j>i ~ 1.45. 

Notice that for dd = md and mc 3> m„, V^J acquires 
the form presented in Eq. (|4.5p . 

But again, dd = mb drives pattern II to a six texture 
zero "ansatz" with 6*13 = 6*23 = 0, implying 

once more a mandatory fine tunning for the free pa- 
rameter dd- 

The random numerical analysis for pattern II throws 
the following numbers (mass parameters in GeV's and 
angles in radians 

as before): 



0.97437 0.22493 0.00250 \ 
0.22478 0.97363 0.03900 , 
0.00856 0.03813 0.99924 J 

numbers again in good agreement with the experimen- 
tal measured values. Finally, the three angles of the uni- 
tary triangle 

of the B decays CP asymmetries are calculated for this 
pattern to be: 

(a, /?, 7)if/^ = (86.78, 16.11, 77.11), 

which again close the triangle and are such that a and 
7 agree with the measured value at la with a 
value for /3 at 2a. 



mu = 0.0023, rUc = 0.560, m* = 172, 



3. Pattern III 



md 0.0029, ms = 0.060, mb = 2.89, 



Proceeding as in the two previous patterns and using 
the notation introduced in appendix (B) we have now 
V^^j^j^j = 0™C/;^"C/J^'^[0;|'^]'^, where the elements of the 
CKM mixing matrix can be expressed as 



dd = 2.8855 < mb, 0i = 1.6, 
which imply the following 3x3 mixing matrix: 



{^CKM)lm 



{or)ii{of)mi + e^'^HOT 

+e"^'{Or)i3iOr)m3 



)l2{Od )vi2 

(4.8) 
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where (f>i = (a" — a™) and 4>2 ~ (a" — + /3" — /3™), which allows us to write the following analytic forms: 



*^ lid 



(mt - rric) 



(mt - mu){mc + m„)(m6 - md)(md - m^) 



dd{mt — mc + ruu) 



i</,2 j mumd{mt + mu){mc — mu){mb — dd)(dd + m^) 
dd{mt — nic + mu){rnt — nic) 



v: 



{mt — rric) 



(mt — m„)(mc + mu)(mi, + ms){md + ms) 



mtmcmtmd{dd + rUs 
dd{mt - mc + m„) 



-y/m^ms{dd + ms)- 



.4,2 I raums{mt + mu){mc - m„)(mj, - dd)(dd + rwd) 
ddirrif — mc + mu)(rnt — nic) 



Vub 



{mt — mc) 



{mt — mu){mc + m„)(m6 + ms){mb + md) 



mtmcmsmd{mt — da) 
dd{mt — mc + m„) 



+ e"''^^/m„m6(mi, - dd) 



4,2 I mumb{mt + Tn„)(mc — mu){dd + ms){dd — md) 
dd{mt — mc + m„)(mt — mc) 



{mt + : 



{mt + mc){mc + m„)(m(, — md){md + m^) 



mtmumbms{dd — ?nd) 
dd{mt ~ mc + mu) 



I 73 T i,*, mcmd{mt - mc){mc - mu){mb - dd){dd + ms) 

\/ mcmd\ad — nrid) — e t — ; r- ; r 

y dd{mt — mc + mu){mt +mu) 



yiii ^ 

* r.s 



{mt + mu) 



{mt + mc){mc + m„)(m(, + ms){md + ma) 



mtmumbmd{dd + ms) 
dd{mt — mc + mu) 



, id,-. I J-, — ; r , i<A2 / mcms{mt - mc){mc - m„)(m(, - dd){dd - ma) 

+ e ^mcms{dd + ms) + e '^^i / — ; -7 ; r 

y dd{mt — mc + mu){mt + mu) 



V, 



{mt + mu) 



{mt + mc){mc + 7nu) {mb + ms ) {mb — md) 
\/ mcmb{mb — dd) 



mtmumsmd{mb — dd) 
dd{mt — mc + mu) 



mcmb{mt - mc){mc - mu){dd - md){dd + m^) 



dd{mt — mc + mu){mt + m„) 



Vti 



{mc — mu) 



{mt + mc){mt — mu){ms + md){mb — md) 



mcmumbms{dd — md) 
dd{mt — mc + mu) 



+ e 



I 7- r iA,-, mtmdimt — mc){mt + mu){mb — dd){dd + ms) 

\/mtmd{dd ~ md) - e a ^ — — ; — - 

y dd(mt — mc + mu){mc — mu) 



V/s" = 



(mc — mu) 



{mt + mc){mt — mu)(m(, + ms){ms + md) 



mcmumbmd{dd + ms) 
dd{mt — mc + mu) 



e"*"^ ^/mtms{dd + m^) + 1 



i4,2 jmtms{mt - mc){mt + m„)(mb ~ dd){dd - md) 



dd{mt — mc + mu){mc — mu) 



ytb — 



{mc — mu) 



{mt + mc){mt — mu){mb + ms){mb — md) 



mcmumdms{mb — dd) 
dd{mt — mc + mu) 



+ e"*"^ ^/mtmb{mb — dd) + 1 



14,2 I mtmb{mt — m^{mt + mu){dd + ms){dd — md) 



dd{mt — mc + m„)(mc — m„) 



r 



(4.9) 



Notice that for dd — md and mt ^ mc ^ m„, V^^ acquires the form presented in Eq. (|4.5 
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But again, dd = nib drives pattern III to a six texture 
zero "ansatz". 

Written in the previous form, all the entries for 
includes two different phases, (jfi and (j)2, 

none of them absorbed immediately by a redefinition 
of a single quark field; but since it is a well known fact 
that the 

SM mixing matrix can be parametrized with only one 
single phase, our analysis makes sense only for the fol- 
lowing three 

different cases, which must be studied separately: 

• Case 1: (f)i 7^ 0, (f>2 ^ 0. 

• Case 2: ^0, 0- 

• Case 3: = 02 7^ 0. 

The numerical analysis for pattern III and for = (/)2 
throws the following numbers (mass parameters in GeV's 
and 

angles in radians as before): 



0.0023, mc = 0.560, mt = 172, 



0.97437 0.22493 0.00247 \ 
0.22408 0.97365 0.03837 , 
0.00817 0.03757 0.99926 / 

numbers again in good agreement with the experimen- 
tal measured values. Finally, the three angles of the uni- 
tary triangle of the B decays CP asymmetries calculated 
for this Pattern are: 



(«: /3> 7)i"'' = (92.73, 16.22, 71.06), 

which again close the triangle and are such that a and 
7 agree with the measured value at \a and /3 at 2a. 



4. Pattern IV 

Proceeding in a similar way we have for this pat- 
tern that the CKM mixing matrix reads now Vq^j^j — 
OuU^''U^"[0''/f, and the elements of the CKM mixing 
matrix can be expressed as 



rud = 0.0029, rus = 0.06, nib = 2.89, 

dd = 2.889, 01 - 02 = 1.6, 
which imply the following 3x3 mixing matrix: 



(4.10) 



where 0i = (a^-a^) and 02 = {ai + K^a^), which 
allows us to write the following 
analytic forms: 
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^iv 



V... 



^mtiricms {mt — dn + rric) + e'"'*'^ •\/mdm„(d„ — nic + m„)(mt + rric — du){irit — rUu — du) 
%J {nid + ms)(~du + TTic + mt){mc + mu)(—du + rric + nit — mu)(mt + rriu) 



— y/mtmcmd {nit — du + nic) + e"^^ -J msmu{du — rric + ?n„)(mt + rric — du}{mt — m„ — du) 
(rrid + ms)(— d„ + + mt){mc + mu){—du + nic + rrit — my){mt + m^) 



,riv _ i02 {mt-du + mc) 



{rric + niu) {mt + m„) 



IV _ —e"^^mc^md{mt — du + mc){du — rric + rriu) + \/msm.cm.umt{mt — du ~ rriu) 

V cd — ^ ; 

^mc{ind + ms){mc + mu){mt — mc){mt + nic — rUu — du) 



yiv _ e"*'^mc\/ms{mt — du + nic){du — rUc + niu) + \/rndnicmumt{mt — du — niu) (4 11) 

y^mc{md + nis){mc + niu){mt — mc){mt + rUc — rUu — du) 
.riv _ j4,-2 \/mc (mt - du +mu) 

^ch 



\J {mt — nic) {rric + niu) 
^/ mdmt{mt — du + mc){du — mc + mu){mt — du — m^) + y^rn^m^^rnZ{du — mc + mu) 
\J {md + nis){du — mc + mu){mt + mc — mu — du){mt — mc){nit + mu) 



e"^'^ msmt{mt — du + mc){du — mc + mu){mt — du — mu) + yjrni^ni^m^{~du + mc — mu) 
\/ {md + ms){du — mc + mu){mt +mc~ mu — du){mt — mc){mt + mu) 
\J mt {du — mc + mu) 
{mt — mc) {mt + mu) ' 



where again the common phase e*"^^ in the last column 
can be removed by a redefinition of the b quark field. 

The numerical analysis for pattern IV arc: (mass pa- 
rameters in GeV's and angles in radians as before): 

m„ = 0.00127, nic = 0.619, nit = 171.7, 



nid = 0.0028, = 0.055, nib = 2.89, 



du = 171.39, 01 = 1.9, 
which imply the following 3x3 mixing matrix: 

0.97331 0.22549 0.00333 \ 
0.22934 0.97242 0.04144 , 
0.00889 0.04163 0.99909 / 

numbers again in good agreement with the experimen- 
tal measured values. 



Finally, the three angles of the unitary triangle of the 
B decays CP asymmetries calculated for this Pattern are: 



(a, /?, 7)^^"^' = (95.23, 19.36, 65.42), 

which not only close the triangle, but the three agree 
with the measured value at Icr. 

5. Pattern V 

Proceeding in a similar way we have for this pattern 
that the CKM mixing matrix reads now 

V^j^j^j = Oi'^U^'^U'^^'^p'^Y , and the elements of the 
CKM mixing matrix can be expressed as 



{V^KM)im = {Ot%,{OfU,+e^^^{Ot%2{OT),n2 
+e'*HOf)n(Of)™3, (4.12) 

where (pi = — ajf and (f)2 = (a.f — a^" — /3J), which 
allows us to write the following analytic forms: 
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' us 



uh 



V 



a/ rribmsmt {dd - md) - e"^^ ^ (^^^ _ ^f^) (^^^ _|_ ,7^^) mdTO„ 
^drf (mfc - m^) (m^ + m^) (»7it + m„) 



y _ yjrtibmdmt {dd + iris) + e"^^ ^ {irib - dd) {dd - md) msrUu 

* 11 



yj {mt + rriu) dd {rm + m^) (m^ + md) 



yjmdmsmt (mf, - dd) + e"^^ yjmbmu{ms + dd) {dd - md) 
a/ (mt + m„) (mfc - m^) (mb + vis) 



\/md{dd 



md) 



y/{mb - md){md + m^) ' 



{mb + ms){md + ms) 



yjmb{mb - dd) 
^(mb - md){mb + m^) 





rtibms {dd - 


' 77id) - e**2 a/ mtmd{mb 


- dd) (dd 


+ ms) 






mu) dd {mb - md) (ms H 


h ?71d) 






mtmd {dd -f 


- ms) + e*'*2 y^mtms{mb 


- dd) (dd 


- md) 






mu) dd {mb + ms) {ms H 


- "Id) 




-\Jm^ 




-dd)+e"l>^^mtmb{dd^ 


1- m.s) {dd 


- md) 






mu) dd {mb + ms) {mb ~ 


- md) 





where the phase 01 can be removed by a redefinition 
of the charm quark field. 

The numerical analysis throws now the following num- 
bers (mass parameters in GeV's and angles in radians as 
before) : 



mu = 0.00127, mc = 0.619, m* = 171.7, 

md = 0.0029, ms = 0.055, mb = 2.89, 

dd = 2.885, (f)2 = 1.9, 
which imply the following 3x3 mixing matrix: 



0.97463 0.22377 0.00266 \ 
0.22361 0.97381 0.04122 , 
0.00891 0.04034 0.99915 J 

numbers again in quite good agreement with the ex- 
perimental measured values. Finally, the three angles of 
the unitary 

triangle of the B decays CP asymmetries calculated for 
this pattern are: 



(a, /3, 7)1^^ = (88.83, 16.36, 74.81), 



which close the triangle and are such that a and 7 
agree with the measured value at Icr, with a value 
for /? at 2a. 



6. Related patterns 

The analysis shows that a weak basis transformation 
using the unitary matrix 



1 
J7 = I 1 I = ^7^ 
1 

shows that UM'^U = A-r^ and also that UM^fU 
M'^'^. So, a pattern like 



M, 



VI 




M, 



VI 



ad 
dd bd 



Cd 



should be physically equivalent to Pattern II. 
A pattern like 



VII 




MY" = 
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should be physically equivalent to Pattern IV. 
And a pattern like 



VIII 



a„ 
al bu 
Cu 



M, 



VIII 



ad 
Cd bd 
a*, b*, dd 



should be physically equivalent to Pattern V. 



V. DISCUSSION OF THE RESULTS 

All the five texture zeros patterns studied here, pro- 
duce numerical results for the CKM mixing matrix and 
CP violating parameters (a, /3, 7) in agreement with the 
measured values up to 2a. Pattern IV reproduces them 
up to 1(7. Since the results are presented for the best fit 
values, all the patterns can be taken predictive, as far as 
the up to date measured numbers are concerned. 

But: how sensitive is our numerical analysis to the 
variations of the quark mass values? Since the answer to 
this question is pattern dependent, let us do it for pattern 
IV, the most predictive one. 

Since mj 3> nic S> m„ and the experimental uncer- 
tainty in the top quark mass Amt is such that Am^ > 
TOc > iriu', then for the up quark sector it is enough to 
consider variations to the top quark mass in the interval 
168.7 < mt{GeV) < 174.7. 

Our results show that the numerical predicted values 
remain stable for the entire interval, as far as the dif- 
ference {rrit — du) stays constant, in a range from 0.28 
to 0.33 [notice that all the CKM entries for pattern IV 
but V/ij^ ^ 1 in Eq. (|4.1ip . depend upon the difference 
(mt - du)]. 

Now, for the quark masses in the down sector, the 
CKM mixing matrix for pattern IV does not depend upon 
the b quark mass (due to the block form of the down 
quark mass matrix in this particular pattern). For the 
other down quark masses we have that the numerical 
results are stable as far as the strange quark mass is in 
the range 0.051 < ms(GeV) < 0.060 and for the down 
quark mass 0.025 < mrf(GeV) < 0.030. 

Can we relate our different parametrization of the 
CKM mixing matrix with other different parameteriza- 
tion? In special to the advocated standard parametriza- 
tion in matrix (j2.6p ?. the answer to this question is an 
academic exercise without much physical sense, because 
the adoption of a parametrization of the flavor mixing is 
an arbitrary fact and not a physical issue. In particular, 
our phase 0i in pattern IV and the phase S of the stan- 
dard parametrization, are not physical quantities. The 
real physical parameters are in this case the CP violating 
mixing angles (a, /?, 7). 

Also, the rotation parameters 9i2, O23 and ^13 of the 
standard parameterization, are not physical entities ei- 
ther. 

For example, in the standard parameterization, 612 = 
t'An~^{Vus/Vud) which is a real number. In any of our 



parametrization, (Vus/Vud) is a complex number, with- 
out direct relation to 9i2 . 

Since in general grounds one expects that the flavor 
mixing parameters do depend on the elements of the 
quark mass matrices, we feel that our CKM mixing en- 
tries arc more realistic that any ad-hoc parametrization. 

Can we do the mathematical exercise to relate the 
phase (jji in pattern IV to the phase S of the standard 
parametrization? (or to any other phase of an arbitrary 
parametrization?). The answer is yes and the simplest 
way to do it is by the use of the Jarlskog [22| invariant, 
which by the sake of convenience we use as: 



j=±im.(K.K6K;K:), 



(5.1) 



which for the standard parametrization takes the form: 
J = ci2C?^C23Si2 513523 sin(5, with a numerical evaluation 
producing [l3l 

(2.80 < J < 3.16) X 10-5. 

For pattern IV, an analytic expression for J can be 
obtained, but it is prohibitive long to quote it here. A 
numerical evaluation is more viable, and for the best fit 
values for the quark masses we get J = 2.87 x 10"^ which 
falls in the allowed range of the standard parametriza- 
tion. Of course, this value is sensitive to the variation 
of the quark mass values. In particular we see that an 
strange quark mass (keeping the other masses in the best 
fit values) in the interval 0.048 < ms(GeV) < 0.056 is 
compatible with the range of the Jarlskog invariant of 
the standard parametrization. 



VI. CONCLUSIONS 

In this paper we have presented our exhaustive study 
of all the possible five texture zeros quark mass matrices 
and their 

connection with the CP violation phenomena for three 
SM families. Analytic and numeric studies were per- 
formed in full detail. 

Let us quote some concluding remarks: 

1. In the context of the SM or its extensions without 
right-handed currents, the quark mass matrices 

Mu and Md can be taken to be Hermitian without loss 
of generality. Non Hermitian quark mass matrices are 
relevant only 

when physics beyond the SM is being considered. 

2. By counting free parameters we have that three 
texture zeros or less in the hermitian quark mass matrices 
of the SM, do 

not imply physical predictions for the elements of the 
flavor mixing matrix, because these texture zeros can al- 
ways be obtained 

in a trivial way by using weak basis transformations. 

3. Each one of the five RRR patterns studied here 
includes two physical relationships between the quark 
masses and the mixing 
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angles, as can be seen from V^^^^; # = /, //, . . . ,V. 
These predictions are case dependent, but far from being 
trivial. 

4. Four and six texture zeros imply one and three 
physical relationships respectively, and allow us to write 
one or the 

three mixing angles as functions of the six quarks 
masses and the CP violation phase. But six texture zeros 
are already 

ruled out by analytic and numerical considerations [Tsl . 

m- 

5. More than six texture zeros are not possible because 
they imply either DetAf^ = or a degenerate quark mass 
spectrum, 

both situations incompatible with the real world. 

6. For the five patterns studied here, number IV (or 
equivalent number VII) reproduce all the experimental 
measured 

values, including the CP violation asymmetries, at la. 
For the other patterns, the /3 angle appears 

at 2a. So, five texture zeros are far from being ruled 
out. 

7. In our analysis for pattern IV, the several con- 
straints for the strange quark mass imply 0.051 < 
ms(GeV) < 0.056. 
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Appendix A: Mass values 

For the experimental mass values used for carrying 
out the numerical calculations in the main text, we have 
adopted the 

following ranges of quark masses [1^ at the Mz energy 
scale (where the Vckm matrix elements are measured). 



Up sector 


Down sector 


mt = 171.7 ±3.0 GeV 
nic = 0.619 ± 0.084 GcV 
m„ = 1.27t^-^2 MeV 


?7i6 2.89 ± 0.009 GeV 

rus = MeV 
md = 2.90ti i9 MeV 



Appendix B: Orthogonal Matrices 

In this appendix we derive the analytic orthogonal ma- 
trices which diagonalize the several 3x3 real quark mass 
matrices coming from the five different RRR patterns. As 
can be seen, the ten quark mass matrices in Table pil[) 
can be grouped in only five different forms. 

For the analysis we write the parameters of the sev- 
eral real mass matrices as functions of the quark mass 
eigenvalues mi, m2 and TO3, with the hierarchy \mi\ < 
\m2\ < {msl, by making use of the following invariant 
forms: tT[Mq], tr[(Afq)2], and det[M5]. 



1. Bi-diagonal form 

For this form the hcrmitian quark mass matrix is 



M, 



dd 




(Bl) 



where the non-zero elements are only in the two diag- 
onals. From the Table, only matrix belongs to this 
form. 

Since det[M^''] = -Cq\aq\^ < for Cq > 0, we use for 
the diagonal quark mass matrix Dg(mi, 7712, —m^). 

Now using Uq'^ = Dg.{\, 1, e'"') we have 



Ijddj^jddpa 







(B2) 



The 3 matrix invariants provide us with the following 
set of equations 



dq -1712+ "I3, 



qi = mim2m3, 

= mim2 -I- 77127713 — 77117713, 



(B3) 



Cqdq 



which allow us to write 



7771 - "^3, 



7772, \aq\ = 17773 (B4) 



with the exact diagonalizing transformation 0^"^^ given 



by 



The light quark masses 77i„ , 771^ and 771^ can be further 
constrained using the mass ratios (23j 



m„/777d = 0.553 ± 0.043 77is/777d = 18.9 ± 0.8 . (Al) 
Notice also that due to the experimental errors 



77l4 ± 77lc ± Till 



771*. 



o 



dd 







± 

1 









(B5) 



which implies a mixing between the first and third 
quark families proportional to 



sm fyi3 = 



m.i 



T7i -|- 7773 
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2. Block Form 

For this form, the hermitian quark niass matrix is: 



q 







flo e 







la^le'"^ 

Cq 

d„ 



(B6) 



From the Table, matrices M^, M^^ and belong 
to this form. 

Using U^f^f = i:>5.(l, e*"^ 1) we have 



|a,| 



'91 '-9 

d, 







(B7) 



We use for the diagonal quark mass matrix 
Dg(mi, —77X2, 7B3). The 3x3 matrix invariants combined 
with the hierarchy \mi\ < \m2\ < jmal allow us to write 



nil - m2, 



(B8) 



\aq\ ^ y/mim2, 



with the solution associated with the exact diagonal- 



izing transformation O^q'^'^ given now by 
/ 



V 



mi4-m2 
mi 
?7ii+m2 





± 



mi 4-"i2 

m2 
m.i+m2 







1/ 



(B9) 



3. Nearest Neighbor form 

For this form, the hermitian quark mass matrix is: 



|6,|e-< 



(BIO) 



From the Tabic, only M^^^ belong to this form. 
Using t/"" = Dg.{l,e< ,e'''<+^^^) we have 



Using the definitions 




A/rnn 



(Bll) 



M' 



q{nn) 



[O7Y M;(fl)07 = Diag(mi, -m2, mg), 

(B12) 

where the sub-indices 1,2,3 in the diagonal forms refer 
respectively to the masses for the quarks u, c and t for 
the up sector, as well as d, s and b for the down sector. 

The three matrix invariants allow us to write now: 



dq = mi - m2 + ms, 
2 _ TO1TO2TO3 

TOl — 7712 + "7.3 



(B13) 



\bq 



2 _ (7(73 - 7772)(7773 + 777l)(7772 - 7^1) 
7771 — 77l2 + 77I3 



and the exact diagonalizing transformation O™^ for 
this particular form is expressed as 



± 



(ms — mi)(m2+mi)(mi — m2+m3) 

(m2+mi ) (7713 +771-2) (mi — m2+m3) 

7niT7i2(m2— TTi-i) 

(m2 +777,3 )(m3 — 777i)(mi —7772+7773) 



7771(7773-7712) 

( 777 3 — 777 1 ) { 777 1 + 771 2 ) 

'^2('^1+T7i3) 

( m 2 + 777 3 ) ( 777 1 + 777 2 ) 

"^3 (7772— 777l) 
( 777 2 + 777 3 ) ( 777 3 — 777 1 ) 

I 





"11 ("12- 


"ll)("ll- 


"I3) 




-"ll){"ll 


^m2)("ii 


-"12 +"13) 




"12 ("13- 


"12){"12- 


"11 ) 


("12 


f mi)(m3 


4-m2)("li 


— m2+"13) 




"13 ("13- 


"12)("11- 


-ma) 


("12 


f "13)("13 


-"ll)("ll 


~"12+"13) 



(B14) 



where one has the freedom to choose two equivalent 
possibilities of phases (the up or down signs). 



J 



For the up quark sector, and due to the fact that 7774 3> 
rric ^ rUu, the former matrix (|B14I) can be expanded as 



(±(1 — r77„c/2) ±V"7-i/c(l — mct/2 — rriuc/'i) ^^y/rrhu{l — rriuc + mct/2) 
±V"^iic(l - muc/2) =f(1 - mct/2 - m„c/2) ±^mct{l - rUuc - mct/2) 
iruct^rn^ ±^/rn^{l — m^c/2 — mct/2) ±{1~ mct/2) 



(B15) 
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where rriij = rrii/mj, i < j; i, j = 1,2,3 
respectively, and we have taken m^t = 0. 



The matrix invariants allow now to write: 



4. Weak basis form 

For this form, the hermitian quark mass matrix is: 



M. 



wb 






ls|e^"' 












(B16) 



From the Table, matrices M^, Mj^, A/j^^ and 
belong to this form. 

Using U^'' ^ £15.(1, e*"™,e'K+'3,"')) wc have 



mi — m2 + 7713 

TOim2TO3 



2 _ {dq - mi){dq + m2)(m3 - dg) 

9 ~ ' 



(B18) 



1^ 



^M'^Cj^y (B17) pressed now as \2A 



where, we use for the diagonal quark mass matrix 
Dg(mi, — m2, m3). Notice here that, for dq = 771,3, bq ~ 

and thus M™'' acquires the block form with an extra 
texture zero. 

The exact diagonalizing transformation can be cx- 



^1 



± 



m2m3(iig — mi) 
(ms — mi)(m2+mi )dq 
mi (dg — mi ) 
(m3 — mi)(mi+m2) 
f^l (•^3 —dq){m2+dq) 
(ms— mi)(mi+m2)dq 



=F 

zb 



mi m3 (dq + m2 ) 
(m2+mi )(m3+m2)dq 



'm2{dq+7n2) 
m2+m3)(mi +"1-2) 



/ m2(m3— dg)(dg— mij 
Y (»n2+"ll)("l3+™'2)dg 



± 

± 
± 



mim2(m3— dg) 
( m2+m3)(m3— mi)dg 
/ m3(m3-dg) 
V ('Ti2+"'3)(™3-"'l) 
/ m3(dg— mi)(iig+m2) . 
Y (»n2+"i3)(»"3— "iljdii / 



(B19) 



5. See-saw form 

For this form, the hermitian quark mass matrix is: 



The 3x3 matrix invariants allow us to write 








(B20) 



From the Table, only matrix belongs to this 

form. Now, for the diagonal quark mass matrix we use: 
Dg(-mi,?7i2,m3). 

Using now U^^ = Dg.(l,e^^°'i^/^q\e'^°'i) we have 




M, 



q{R)- 



(B21) 



2 _ (m3 + m2 - dq){dq -7712 + mi){-dq + - mi) 

— dq — TOi + m2 + 7713 

Cq = -dq - mi + ?7l2 + ?7l3, 
2 771177127713 

I"?! " ~i \ ; — ' 

— dq — 7)71 + 7712 + 7713 

(B22) 

where notice that for dq = m^ — mi, bq = and M^^ 
becomes the two diagonal form with an extra texture 
zero. 

The exact diagonalizing transformation O^" is given 
now by 



"^2^3 ( — dq4-m2+m3) 
(mi+m2)(mi +m3)( — dq — m.i+m2+m3) 

mims ( — dq+"T,3— mi ) 
(mi+m2)( — dq— mi+m2 4-m3)( — m2+m3) 

mi m2 (dq — m2 + mi ^ 
Y Y ( — m2-|-m3)( — dq— 7ni+m2+m3)(7?ii-|-m3) 



mi (dq — m2 + mi )( — dg+ms + mi ) 
(7?il+m2)(mi+m3)( — dq — 7ni+m2+m3) 

m2{dq-^mi~m2)[~dq+m2 + m-i^ 
(mi+m2)( — dg — mi+m2+m3)( — m2+m3) 

"^3 ( — "^g +"^2 +"^3 ) ( — — +"T'3 ) 
( — dg— mi+m2+m3)( — m2+m3)(mi + m3) 



/mi( 


-dq+rri2+m3) 


(mi- 


-m2)(mi+m3) 


m2( 


-dq— mi+m3) 



( 7^^. 1 + Tn 2 ) ( — 2 + 3 ) 

(dg— m2+mi^m3 
(-m2+m3)(mi+m3) J 



(B23) 
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